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Abstract. Based on a comparison of dynamical models with observations of the interstellar gas in 6 Virgo cluster 
spiral galaxies a first complete ram pressure stripping time sequence has been established. The observational 
characteristics of the different stages of ram pressure stripping are presented. The dynamical models yield the 3D 
(— I ' velocity vectors of the galaxies, peak ram pressures, and times to peak ram pressure. In the case of a smooth, 

^ static, and spherical intracluster medium, peak ram pressure occurs during the galaxy's closest approach to the 

cluster center, i.e. when the galaxy's velocity vector is perpendicular to its distance vector from the cluster center 
(M87). Assuming this condition the galaxy's present line-of-sight distance and its 3D position during peak ram 
pressure can be calculated. The linear orbital segments derived in this way together with the intracluster medium 
^ ' density distribution derived from X-ray observations give estimates of the ram pressure that are on average 

a factor of 2 higher than derived from the dynamical simulations for NGC 4501, NGC 4330, and NGC 4569. 
Resolving this discrepancy would require either a 2 times higher intracluster medium density than derived from 
X-ray observations, or a 2 times higher stripping efficiency than assumed by the dynamical models. Compared to 
NGC 4501, NGC 4330, and NGC 4569, NGC 4388 requires a stiU 2 times higher local intracluster medium density 
or a direction which is moderately different from that derived from the dynamical model. A possible scenario for 
the dynamical evolution of NGC 4438 and M 86 within the Virgo cluster is presented. 
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1. Introduction tions of these clusters are characterized by the presence 

of a large number of post starburst galaxies. Poggianti 

The evolution of a spiral galaxy is strongly affected by ^^ggg) concluded that the most evident effect due 

q . Its environment givmg rise to the density-morphology re- ^^^^^^^ environment is the quenching of star for- 

_C3 : lation (Dressier 1980): early type galaxies are more fre- ^^^-^^ ^^^^^^ ^j^^^ enhancement. They found two dif- 

quently found m high density environments, i.e. galaxy ^^^^^^^ g^j^^y evolution timescales in clusters: (i) a rapid 

clusters. At intermediate redshifts (z < 0.5) the fraction of ^.^^^^^j ^j^^ ^ ^^^^ ^ ^1^^ transforma- 

blue starformmg galaxies m clusters is higher than that m ^.^^ morphology (several Gyr). A recent study of two 

local clusters (Butcher-Oemler effect; Butcher & Oemler ^^^^j^^^ intermediate-redshift galaxy clusters (Moran et 

1978, 1984). Together with the resuff that the fraction of 2007) confirmed these findings: spiral galaxies within 

SO galaxies increases with decreasing redshift (Dressier et ij^falling groups already begin a slow process of conver- 

al. 1997; Fasano et al. 2000), this indicates that infalling ^-^^ -^^^ gQ^^ q^^^ ^^^^ ^^^^^ ^1^^^^^ ^^^^^ p^^^, 

spiral galaxies are transformed into SO galaxies. Since mor- ^^^^ stripping dominates where the intracluster medium is 

phology IS closely related to the star formation history of a Otherwise, the slow transformation process goes on. 

galaxy, it is not surprising that the average galaxy proper- Whereas the slow transformation is most probably due to 

ties related to star formation also depend on local density gravitational interactions (galaxy-cluster potential, Byrd 

(Hashimoto et al. 1998; Lewis et al. 2002; Gomez et al. ^ Valtonen 1990, Valluri 1993; galaxy-galaxy; galaxy ha- 

2003; Kauffmann et al. 2004; Balogh et al. 2004). The anal- .^ssment, Moore et al. 1996, 1998), the rapid timescale is 

ysis of galaxies in ten distant clusters (Dressier et al. 1999; p^^g^^^^ stripping. 
Poggianti et al. 1999) has shown that the galaxy popula- 

The Virgo cluster of galaxies represents an ideal labo- 

Send offprint requests to: B. Vollmer: bvollmer@astro.u- ratory to study the effects of ram pressure stripping in de- 

strasbg.fr tail. It is the only cluster in the northern hemisphere where 
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one can observe the ISM distribution and kinematics of 
cluster galaxies at a kpc resolution (1 kpc ~ 12"j3- The 
Virgo cluster is dynamically young and spiral-rich. Most of 
the clusters spiral galaxies are Hi deficient, i.e. they have 
lost a significant amount of their ISM (Chamaraux et al. 
1980, Giovanelli & Haynes 1983). Imaging Hi observations 
have shown that these galaxies have truncated Hi disks ( 
GiovaneUi & Haynes 1983, Cayatte et al. 1990). Thus, the 
cluster environment changes the Hi content and morphol- 
ogy of Virgo cluster spiral galaxies. However, Virgo spiral 
galaxies are not CO-deficient (Kenney & Young 1989). 
In a new Hi imaging survey of Virgo galaxies (VIVA: 
VLA Imaging of Virgo galaxies in Atomic gas), Chung 
et al. (2007) found seven spiral galaxies with long Hi tails. 
These galaxies are found in intermediate- to low-density 
regions (0.6-1 Mpc in projection from M87). The tails are 
all pointing roughly away from M87, suggesting that these 
tails are due to ram pressure stripping. Therefore, ram 
pressure stripping already begins to affect spiral galaxies 
around the cluster Virial radius. The further evolution of a 
galaxy depends critically on its orbit (see, e.g., VoUmer et 
al. 2001), i.e. a highly eccentric orbit will lead the galaxy 
at a high velocity into the cluster core, where the intra- 
cluster medium is densest and ram pressure will be very 
strong. 

In this article I will use the results from about 10 
years of studying the effects on ram pressure stripping 
on individual Virgo cluster galaxies to (i) construct a first 
complete model-based time sequence and (ii) to verify if 
the model-deduced galaxy orbital segments are consistent 
with observations. In Sect. 2 the method to derive the 
ram pressure stripping parameters is described and the 
uncertainties in the parameters are given, the ram pres- 
sure stripping time sequence is established, followed by a 
description of the observational characteristics of the dif- 
ferent stages of ram pressure stripping. In Sect. [3] it is 
investigated if the model-derived parameters are consis- 
tent with eccentric galaxy orbits within the Virgo cluster. 
Sect. [His devoted to the ram pressure stripping efhciency, 
and the conclusions are given in Sect. [5] 

2. The first complete model-based ram pressure 
stripping time sequence 

Within the last 6 years a sample of 6 Virgo cluster spi- 
ral galaxies was investigated in detail (Table [1] contains 
their B magnitudes, radial velocities, optical radii, and dis- 
tances to M 87; for references see Table [5]) based on VLA 
Hi and IRAM 30m CO observations, and the dynamical 
models described in VoUmer et al. (2001). These galax- 
ies are located at projected distances of 0.3-1.0 Mpc from 
the cluster center, M 87 (Fig. [3]). Besides tell-tale signs in 
their Hi properties (see Sect. 2.4), all galaxies show asym- 
metric ridges of polarized radio continuum emission in the 
outer galactic disk (VoUmer et al. 2007), which are tracers 
of interactions with their environment and can be used as 

^ I use a distance of 17 Mpc to the Virgo cluster 



Table 1. Sample of Virgo spiral galaxies. 



Name 


B 




R25 






(mag) 


(kms^^) 


(') 


(Mpc) 


NGC 4330 


14.0 


1572 


4.0 


0.6 


NGC 4388 


12.2 


2515 


5.5 


0.4 


NGC 4438 


12.0 


98 


8.9 


0.3 


NGC 4501 


10.6 


2270 


6.8 


0.6 


NGC 4522 


13.6 


2324 


3.5 


1.0 


NGC 4569 


11.8 


643 


10.2 


0.5 



diagnostic tools to determine the kind of interaction (grav- 
itational or ram pressure) and the interaction parameters 
(VoUmer et al. 2006; Soida et al. 2006). 

2.1. Determining the stripping parameters 

In the following I describe how the stripping parameters 
are determined by comparing the dynamical simulations 
with observations. For a direct comparison between model 
and observations of a given galaxy the observed parame- 
ters are: its systemic radial velocity, projected distance 
from the cluster center, disk inclination angle, disk posi- 
tion angle, gas distribution, and velocity field. With the 
help of the dynamical model the maximum ram pressure 
Prani^i the time to maximum ram pressure trps, and the an- 
gle i between the galactic disk and the ram pressure wind 
direction were determined in the following way: 

1. Maximum ram pressure and ram pressure wind incli- 
nation angle: 

the galaxy's Hi deficiency {\og{M^^ / M^^)) , where 
is the expected and Mgj^ the observed Hi mass, 
sets the quantity sini: 

M^V<i « sin\^{^ + 10°)) (1) 

(Eq. 21 of VoUmer et al. 2001), where p™^^ is the max- 
imum ram pressure and M^/^ the initial and final Hi 
masses of the simulation, respectively. Typically, the 
Hi distribution of a stripped spiral galaxy shows a 
sharp outer edge on one side and more extended, ex- 
traplanar gas on the opposite side. The determination 
of the inclination angle i is mainly based on the mor- 
phology of the extraplanar gas. The derived maximum 
ram pressure is generally close to that derived from the 
Gunn & Gott criterion (Gunn & Gott 1972): 

Pram = SiSM^rot^"^ : (2) 

where R is the galactocentric distance, Wrot the rota- 
tion velocity, and Sism ~ 10^^ cm~^ the ISM surface 
density at the sharp edge of the Hi distribution. 

2. Temporal ram pressure profile: 

the profile is fixed by the maximum ram pressure due 
to the galaxy motion within the intracluster medium 
. VoUmer et al. (2001, 2008) calculated temporal ram 
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pressure profiles for a galaxy on eccentric orbits within 
the gravitational potential of the Virgo cluster as- 
suming an intracluster medium density derived from 
ROSAT X-ray observations (Schindler et al. 1999). 
Their correlation between the maximum ram pressure 
and the width of the temporal stripping profiles is used 
for the dynamical models. 

3. Azimuthal viewing angle: 

for the comparison with observations the simulated 
galaxy has to be projected onto the sky. The major 
axis position angle and inclination of the galaxy de- 
fine a plane in three dimensional space. The simulated 
galaxy can then be rotated within this plane by an az- 
imuthal viewing angle (see e.g. VoUmer et al. 2008). 
The three-dimensional model wind direction, the line- 
of-sight velocity of the galaxy, and the projected ICM 
wind direction can thus all be expressed as functions 
of the azimuthal viewing angle. 

4. Time to, or from ram pressure peak: 

finally, the model snapshot is chosen which reproduces 
best the gas morphology, velocity field, and observed 
line-of-sight velocity of the galaxy. This sets the time 
to or from ram pressure maximum, i.e. the stripping 
age. In practice, ram pressure simulations are done for 
several different wind directions. 

2.2. Uncertainties of the model-derived interaction 
parameters 

The uncertainties in the model-derived interaction param- 
eters can only be estimated, since the choice of the "best- 
fit" model is made by eye and only a few simulations 
are made. For models where the galaxy has already ex- 
perienced the ram pressure peak (post-peak models), the 
maximum ram pressure has left its imprint on the gas 
distribution and velocity field. On the other hand, in the 
case of a pre-peak model, the peak ram pressure to be at- 
tained depends on assumptions on the galaxy orbit within 
the cluster. In the following I estimate the uncertainties of 
the interaction parameters based on the case of NGC 4522 
(VoUmer et al. 2006) for post-peak, and on NGC 4501 
(VoUmer et al. 2008) for pre-peak ram pressure stripping: 

— post-peak: the uncertainties in the peak ram pressure 
Pram ; time to peak ram pressure irps, and inclination 
angle i between the disk plane and the ram pressure 
wind direction are: ApJ"^ - 20 Myr, At^ps ^0.1 t^ps, 
and Ai - 10°. 

— pre-peak: the uncertainty in i is Ai ~ 10°. Maximum 
ram pressure can only be determined if a galaxy orbit 
is assumed. The uncertainties in p™m and t^-pg are then 
Ap'^^^ ~ 20 Myr and At,ps ~ 0.1 Up^. 

I thus estimate the total uncertainties on the model- 
derived parameters to be Ap™^ = 20 Myr, Airps = 
0.1 trps, and Ai = 10°. It should be noted that these un- 
certainties do not significantly affect the following results 
and conclusions. 



2.3. Observational gas characteristics 

Based on our Hi and CO observations the Virgo cluster 
spiral galaxies can be divided into the following 4 classes, 
based on their gas morphology. These classes can be ar- 
ranged to form a time sequence for ram pressure stripping 
in the Virgo cluster (Fig.[Tl as will be argued hereafter). 

— (i) gas disk truncated near the optical radius {R25) 
and asymmetric outer gas disk with a tail structure 
(NGC 4501, NGC 4330); 

— (ii) strongly truncated gas disk and extraplanar high 
surface density gas (NGC 4522, NGC 4438); 

— (iii) strongly truncated gas disk and low surface bright- 
ness gas (NGC 4388); 

— (iv) truncated gas disk and perturbed low surface 
brightness arms (NGC 4569). 

2.4. The ram pressure stripping time sequence 

The dynamical models tell us that class (i) galaxies are 
observed before peak ram pressure, class (ii) galaxies are 
close to peak ram pressure, and galaxies of classes (iii) and 
(iv) are observed after peak ram pressure (Fig. [T]). 

Growl & Kenney (2008) independently confirmed the 
model-derived stripping ages (time to peak ram pressure) 
by analyzing the stellar populations of the stripped, gas- 
free outer parts of the disks of NGC 4388, NGC 4522, and 
NGC 4569. 

A more detailed description of the different stripping 
stages is: 

— pre-peak (class (i)): 

The galaxies (NGC 4501 and NGC 4330) are at the 
beginning of a ram pressure stripping event (VoUmer 
et al. 2008; VoUmer et al. in prep.), i.e. they are ap- 
proaching the cluster center. The outer gas disk has 
already been removed and the gas of the inner disk is 
just beginning to be affected by ram pressure. 

— near peak (class (n)): 

NGC 4438 is at present undergoing strong ram pres- 
sure stripping. Despite the strong tidal perturbation, 
ram pressure is the dominant effect on the observed 
gas distribution and kinematics (VoUmer et al. 2005). 
NGC 4522 is also close to peak ram pressure. The gas 
disks of both galaxies are strongly truncated and the 
extraplanar gas has a high column density. 

— post-peak (class (iii) and (iv)): 

after peak ram pressure, i.e when the galaxy leaves the 
cluster core again, the removed ISM is accelerated and 
expands. The extraplanar gas tail becomes larger and 
its gas surface density decreases. This is observed in 
NGC 4388 (Oosterloo & van Gorkom 2005) where the 
dynamical model yields a time to peak ram pressure 
of - 100 - 200 Myr. 

About 200 Myr later the gas tail is still expanding, its 
surface density is decreasing, and the gas is evaporated 
by the hot intracluster medium. The gas tail is thus no 
longer detectable in Hi. Ram pressure is decreasing and 
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observations 



model 



NGC 4501 VoUmer et al. (2008) 
NGC 4330 Chung et al. (2007) 

Abramson et al. (in prep.) 
NGC 4438 VoUmer et al. (2005) 
NGC 4522 Kenney et al. (2004) 

VoUmer et al. (2004b) 
NGC 4388 VoUmer & Huchtmeier (2003) 

Oosterloo & van Gorkom (2005) 
NGC 4569 VoUmer et al. (2004a) 



VoUmer et al. (2008) 
VoUmer et al. (in prep.) 

VoUmer et al. (2005) 
VoUmer et al. (2006) 

VoUmer & Huchtmeier (2003) 

VoUmer et al. (2004a) 



the pushed gas, which has not been accelerated to the 
escape velocity, falls back onto the galaxy and begins 
to resettle there. Thus, the only traces of the past in- 
teraction are kinematically perturbed gas arms located 
at the outer edge of the gas disk due to resetthng gas, 
as observed in NGC 4569. 

2.5. How to recognize ram pressure stripped galaxies 

From the work on the 6 galaxies (Table H]) I have identi- 
fied the following observational characteristics ("smoking 
guns") for the different stages of ram pressure stripping: 

— pre-peak ram pressure stripping: 

in the case of the detection of extraplanar gas, in- 
creasing, moderately strong ram pressure (pram ^ 
5000 cm~'^kms~^) leads to extraplanar gas of mod- 
erate mean surface density (several 10""^^ cm~^) close 
to the galactic disk with a continuous velocity field 
from the disk to the extraplanar region (NGC 4501, 
NGC 4330). 

— near peak ram pressure stripping: 

ongoing strong ram pressure stripping 
(Pram ^ 20000 cm~^kms~^) gives rise to high 
mean surface density (several lO^'' cm~^) extraplanar 
gas close to the galactic disk with a continuous 
velocity field from the disk to the extraplanar region 
(NGC 4522, NGC 4438). 

— post-peak ram pressure stripping: 

for intermediate stripping ages (100-200 Myr after 
peak) gas can be found at large distances from the 
galactic disk. The mean surface density of this extra- 
planar gas is low (a few lO^^ cm-^; NGC 4388). The 
velocity field is continuous from the disk to the extra- 
planar region. 

For higher stripping ages (> 200 Myr) perturbed gas 
arms are visible at the outer edge of the gas disk show- 
ing a discontinuous velocity field from the disk to the 
extraplanar region (NGC 4569). 

At the very end of this sequence (> 300 Myr) stripped 
spiral galaxies show a truncated symmetric and unper- 
turbed gas disk. 

Galaxies with stripping ages < 200 Myr show ridges 
of polarized radio continuum emission at the outer edge 



of the disk opposite to the extraplanar gas. These ridges 
are due to ram pressure compression of the large-scale 
magnetic field. Galaxies which are observed more than 
200 Myr after peak ram pressure can also show enhanced 
polarized radio continuum emission due to shear motions 
of re-accreting or resettling gas (Otmianowska-Mazur & 
Vollmer 2003). 

The time window during which one can still identify 
perturbations due to ram pressure stripping is ^ 500 Myr. 
The orbital period of the orbits shown in Fig. [2] is several 
Gyr. This relatively short time window for the detection of 
ram pressure induced perturbations explains the rareness 
of Hi tails in the VIVA sample (7 out of 50 galaxies; Chung 
et al. 2007). 

3. Orbits of ram pressure stripped galaxies 

The ram pressure stripping time sequence presented in 
Sect. 12.41 does not depend on any assumption on galaxy 
orbits. It only relies on the model-derived time to peak 
ram pressure (pre-peak, near-peak, or post-peak). In the 
following, I investigate if the model-derived parameters 
are consistent with realistic galaxy orbits within the Virgo 
cluster, given the positions, velocities and directions of the 
galaxies within the cluster. 

The dynamical model yields the 3D velocity vector of 
the galaxy, the peak ram pressure, and the time to peak 
ram pressure. I assume linear orbits over an interval of 
±200 Myr around peak ram pressure. Fig.[2]shows typical 
orbits in the gravitational potential of the Virgo cluster 
(Schindler et al. 1999), confirming this assumption as a 
zero-order approximation. To calculate distances, I further 
estimate a mean galaxy velocity v on this linear trajectory 
with respect to the cluster mean, which is based on the 
peak ram pressure and the time to peak ram pressure from 
the galaxy orbit simulations of Vollmer et al. (2001) and 
on the simulation shown in Fig.[2l The uncertainty in this 
mean velocity is about Aw ^ 0.2u. 

In the case of a smooth, static, and spherical intraclus- 
ter medium, peak ram pressure occurs at the galaxy's clos- 
est approach to the cluster center, i.e. when the galaxy's 
velocity vector is perpendicular to its distance vector. This 
condition allows the calculation of the galaxy's present 
line-of-sight distance and its 3D position during peak ram 
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NGC 4438 





model 
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NGC 4522 
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NGC 4388 
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"ISOMyr alter peak -SOOMyr after peak 



Fig. 1. Model-based complete ram pressure stripping time sequence for Virgo cluster spiral galaxies. NGC 4501 
(greyscale: Hi, contour: polarized radio continuum emission; VoUmer et al. 2008) and NGC 4330 (greyscale and 
contours: Hi; Chung et al. 2007; VoUmer et al. in prep.) are approaching the cluster center and are thus in a stage 
of pre-peak ram pressure (class (i)). NGC 4522 (greyscale: Hi, contour: polarized radio continuum emission; Kenney 
et al. 2004; VoUmer et al. 2006) and NGC 4438 (observations: CO spectra on optical image; model: greyscale: gas 
surface density; contours: stellar distribution; VoUmer et al. 2005) are close to peak ram pressure (class(ii)). NGC 4388 
(class(iii); observations: contour: Hi Oosterloo & van Gorkom 2005; model: greyscale: gas surface density; contour: 
stellar distribution; VoUmer & Huchtmeier 2003) and NGC 4569 (class (iv); greyscale and contours: Hi; VoUmer et al. 
2004a) are leaving the cluster center. 



pressure (Table [3]0: let the galaxy's 3D velocity vector be 
V and the present position x = {x,y,z), where the line- 
of-sight component z is unknown. Within a time At the 
galaxy moves a distance of As = vAt. The galaxy's po- 
sition at peak ram pressure is then Xq = {xo,yo,zo) — 
X + As. During peak ram pressure the velocity and the 
distance vectors are perpendicular, i.e. f • a;o = 0. In this 
way Zo and subsequently z and the minimum distance to 
the cluster center at peak ram pressure can be calculated. 
Since I estimate the uncertainty in the direction of the 
model galaxy's 3D velocity vector to be ~ 10°, I searched 
for the direction leading to the minimum impact param- 
eter (closest distance to the cluster center) within a 10°- 
cone around the model velocity vector. 

Changing the galaxy's mean velocity by 20% leads to a 
20% variation in the galaxy's 3D position. A change in the 
time to peak ram pressure has the same effect. The largest 
variation of the resulting orbital parameters is observed 
when the galaxy's 3D velocity vector is modified. This 
variation depends strongly on the initial, model-derived 



^ Negative numbers indicate directions to the left, bottom, 
and in front of M 87. 



velocity vector. Whereas the variation in the galaxy's im- 
pact parameter Dmin due to a change of 10° in the veloc- 
ity vector is modest ( < 20%) for NGC 4501, NGC 4438, 
NGC 4388, and NGC 4569, it is a factor of 2 for NGC 4330. 
The resulting galaxy distances from the cluster center 
along the line-of-sight vary by ~ 100 kpc for these four 
galaxies, and by 1 Mpc for NGC 4330. Altogether, I es- 
timate the uncertainty of the galaxy's smallest distance 
to M 87 to be ~ 30%. The linear orbital segments deter- 
mined in this way together with the minimum distances 
to the cluster center are presented in Fig. [3] and [4] 

The galaxies have random velocity vectors within the 
cluster. The range of the line-of-sight positions is some- 
what smaller, but comparable to the ranges along the right 
ascension and declination axes. NGC 4569 has the small- 
est (210 kpc), NGC 4330 the largest (620 kpc) impact 
parameter. 

4. The efficiency of ram pressure stripping 

It is now possible to check the consistency between the or- 
bital segments deduced in Sect. [3] and the maximum ram 
pressure determined from the dynamical models (Table [2] 
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Table 3. Orbital parameters with respect to M 87. Negative numbers indicate directions to the left, bottom, and in 
front of M 87. Row (l)-(3): present galaxy position; row (4)-(6): 3D unit vector of the galaxy's motion within the 
cluster from the dynamical models; row (7)-(9): 3D unit velocity vector leading to the minimum impact parameter 
-Dmin; row (10): total galaxy velocity; row (11): normalized ram pressure Pnorm = (piCMi^|aiaxy)/(Poi'o)i where po—10~'^ 
cm~^ and wo=1000 kms~^; row (12): time to peak ram pressure (negative values are before peak ram pressure); row 
(13)~(15): galaxy position at peak ram pressure assuming a spherical, smooth, and static intracluster medium; row 
(16): impact parameter or distance at peak ram pressure. In Appendix A it is argued that NGC 4438 does most 
probably not follow the given orbit but is rather located ^ 1 Mpc behind M 87. Since NGC 4522 is being stripped by 
a moving intracluster medium, its impact parameter cannot be determined. The uncertainties in the parameters are 
discussed in Sect. 12.21 





NGC 4501 


NGC 4330 


NGC 4438 


NGC 4388 


NGC 4569 


NGC 4522 


X (kpc) 


-86 


559 


227 


374 


-464 


-210 


y (kpc) 


602 


-303 


183 


80 


229 


-955 


2 (kpc) 


-1 


125 


(-308) 


320 


-17 


(500) 




0.48 


0.05 


-0.81 


0.09 


-0.64 


-0.73 




-0.44 


0.99 


-0.20 


-0.75 


0.53 


-0.60 




0.76 


0.08 


-0.55 


0.66 


-0.56 


0.33 




0.44 


0.18 


(-0.72) 


0.24 


-0.75 




Uy 


-0.59 


0.98 


(-0.13) 


-0.65 


0.52 




riz 


0.67 


0.12 


(-0.68) 


0.73 


-0.41 




V (kms^-"") 


1500 


1700 


(2000) 


1700 


1500 


2000 


Pnorm 


20 


20 


50 


50 


40 


50 


A t (Myr) 


-250 


-100 


10 


150 


300 


50 


xo (kpc) 


88 


592 


(242) 


309 


-108 




yo (kpc) 


368 


-128 


(186) 


255 


-15 




zo (kpc) 


266 


146 


(-294) 


127 


178 




Dmin (kpc) 


460 


620 


(420) 


420 


210 





and [3]). The link between these two quantities is the den- 
sity distribution of the intracluster medium, again assum- 
ing that it is smooth, static, and spherical. Schindler et al. 
(1999) determined the intracluster medium density distri- 
bution in the Virgo cluster from ROSAT X-ray observa- 
tions, assuming hydrostatic equilibrium: 

p=pc{i+iy'' (3) 

where tq is the core radius and pc is the central density. 
I assume pc = 4.2 x 10"^ cm^^, rc = 13.4 Mpc, and 
(} = (Schindler et al. 1999). For the galaxy velocity at 
peak ram pressure I assume v = 2000 kms^^ based on the 
galaxy orbit simulations of VoUmer et al. (2001) and the 
simulation shown in Fig. [51 I define the ram pressure effi- 
ciency as a constant ^ with Pram — Cpicmw^, where picM 
is the intracluster medium density and v the galaxy veloc- 
ity with respect to the intracluster medium. In previous 
works ^ = 1 has been assumed. 

The ram pressure Pram — Picmv^ can be seen together 
with the values derived for our 6 sample galaxies in Fig. [5l 
Except for NGC 4569 the assumed intracluster medium 
density distribution leads to a too small maximum ram 
pressure compared to the values derived from the dynam- 
ical models. Enhancing the intracluster medium density 
distribution by a factor of 2 or decreasing /3 = | leads 
to more consistent results. Thus, I conclude that either (i) 



-? 




-3 
-3 



(Mpc) 



Fig. 2. Galaxy orbits in the Virgo cluster with peri- to 
apocenter ratios of about 1:10, giving rise to ram pressure 
maxima between 10 po and 50 pq, where the normalisation 
is Po = pqVq with po=10~'* cm~^ and t>o=1000 kms^^. 
The inner box delimitates the extent of Fig. [31 
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I n4501 



Fig. 3. Orbital segments from the location of peak ram 
pressure to the present location of the 6 Virgo cluster spi- 
ral galaxies. A spherical, smooth, and static intracluster 
medium is assumed. The arrows indicate the galaxies mo- 
tions projected on the sky. Linear orbital segments are 
assumed. The dotted lines connect the location of peak 
ram pressure and the cluster center (M 87). The 3D or- 
bital segments are shown in Fig. [21 



the stripping efficiency and/or (ii) the intracluster medium 
density has been underestimated. 

Compared to NGC 4330, NGC 4501, and NGC 4569, 
the two spirals which are located close to M 86, NGC 4388 
and NGC 4438, require a still 2 times higher peak ram 
pressure than expected from a smooth and static intraclus- 
ter medium assuming a higher stripping efficiency and/or 
a higher intracluster medium density (see Appendix A). 

For NGC 4388, increasing the time from peak ram 
pressure from 150 Myr to 200 Myr as advocated by Growl 
& Kenney (2008) based on an optical spectrum and UV 
flux of the gas-free outer disk, would increase its im- 
pact parameter from 420 kpc to 460 kpc. Decreasing the 
time from peak ram pressure to 100 Myr (VoUmer & 
Huchtmeier 2003), would decrease the impact parameter 
by less than 10%. Allowing for a larger uncertainty in 
the galaxy's velocity vector yields an impact parameter 
of 310 kpc for the direction n = (0.54, -0.48,0.70). The 
angle between this direction and that of Table [3] is 18° in 
the plane of the sky. I conclude that for NGC 4388 either 
the intracluster medium is locally a factor of 2 denser than 
predicted by a spherical, smooth, and static intracluster 
medium distribution, or that the galaxy's unit velocity 
vector is ^ 20° different from the values given in Table [3] 

The cause for the unusual extent of NGC 4388's Hi 
tail and the complicated case of NGC 4438 are discussed in 
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Fig. 4. Three dimensional linear orbital segments of the 
6 Virgo cluster spiral galaxies. The dotted lines connect 
the location of peak ram pressure and the cluster cen- 
ter (M 87). A spherical, smooth, and static intracluster 
medium is assumed. The orbital segment of NGC 4438 is 
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Fig. 5. Normalized 
{piCMWgaiaxy)/(PoWo): ^here po=10 
kms~^) as a function of cluster radius. The solid line 
traces the normalized ram pressure assuming the in- 
tracluster medium density distribution picM (Eq. [3]) 
of Schindler et al. (1999) and v = 2000 kms^^ The 
normalized ram pressure with 2 picM (dashed line) and a 
different (3 (dotted line) are also shown. The distance of 
NGC 4654 is a lower hmit. 



Appendix A. There I advocate a scenario where NGC 4438 
has been stripped by the hot X-ray gas halo associated 
with M 86, has recently left M 86's Mach cone, and is now 
being stripped by the intracluster medium of the Virgo 
cluster. This scenario implies that NGC 4438's orbit is 
different from that given by Tabled Its trajectory is close 
to M 86, i.e. ~ 1 Mpc behind M 87. In this case, NGC 4438 
cannot be used to test the consistency between galaxy 
orbits, dynamical models, and observations of the galaxy's 
ISM. 

For comparison I have also added NGC 4522 and 
NGC 4654 to Fig. 5. As demonstrated in Kenney et al. 
(2004), the peak ram pressure for NGC 4522 is far too 
high to be explained in our smooth and static intracluster 
medium scenario. Indeed, Kenney et al. (2004), VoUmer et 
al. (2004), and VoUmer et al. (2006) claim that NGC 4522 
is encountering a moving intracluster medium which is 
presumably associated to M 49. Therefore, its line-of-sight 
distance cannot be determined. Since the galaxy subclus- 
ter associated with M 49 is falling into the Virgo cluster 
from behind (Irwin & Sarazin 1996; Biller et al. 2004) and 
the intracluster medium of M 49 is interacting with the 
intracluster medium of M87 (Shibata et al. 2001), I have 
assumed a line-of-sight distance of z = 500 Mpc. Based on 
the dynamical model, NGC 4654 experiences a very low 
ongoing ram pressure together with a gravitational inter- 
action which occurred in the recent past (VoUmer 2003). 
Again, its line-of-sight distance cannot be determined and 
the plotted value is a lower limit. The low ram pressure in- 
dicated that NGC 4654 might be located as far as 2 Mpc 
behind the cluster center (M 87). 



I conclude that the efhciency of ram pressure strip- 
ping derived from a sample of 4 Virgo spiral galaxies 
(NGC 4501, NGC 4330, NGC 4388, and NGC 4569) is 
about unity. The ram pressure stripping efficiencies of 
NGC 4522 and NGC 4654 are consistent with being close 
to unity. 



+n4522 5_ Conclusions 



A first complete ram pressure stripping time sequence 
could be established by combining the results of detailed 
comparisons between dynamical models and observations 
of the interstellar medium in ram pressure-stripped galax- 
ies in the Virgo cluster. It is possible to observe ram pres- 
sure induced perturbations ^ 300 Myr around the galaxy's 
closest approach to the cluster center, i.e. when peak ram 
pressure occurs if a spherical, smooth, and static intraclus- 
ter medium distribution is assumed. The relative brevity 
of this period compared to the galaxy's orbital timescale 
(several Gyr) explains the rareness (r^l4%) of observed 
ram pressure induced perturbations in Virgo spiral galax- 
ies. 

Observationally the different stages of ram pressure 
stripping can be recognized in the following way: 

^ increasing, moderately strong ram pressure (> 50 Myr 
before peak): 

moderately truncated Hi disk, extraplanar gas of mod- 
erate surface density, continuous velocity field between 
the disk and the extraplanar region, ridge of polarized 
radio continuum emission at the outer gas disk oppo- 
site to the extraplanar region. 

— ongoing strong ram pressure (near peak): 

strongly truncated Hi disk, extraplanar gas of high sur- 
face density, continuous velocity field between the disk 
and the extraplanar region, ridge of polarized radio 
continuum emission at the outer gas disk opposite to 
the extraplanar region. 

— decreasing ram pressure (< 200 Myr after peak): 
strongly truncated Hi disk, extended extraplanar gas 
of low surface density, continuous velocity field be- 
tween the disk and the extraplanar region, ridge of 
polarized radio continuum emission at the outer gas 
disk opposite to the extraplanar region. 

— decreasing ram pressure (> 200 Myr after peak): 
strongly truncated Hi disk, perturbed outer gas arms, 
discontinuous velocity field between the disk and the 
extraplanar region, ridge of polarized radio continuum 
emission at the outer gas disk opposite to the extrapla- 
nar region, possible ridge of polarized radio continuum 
emission at the outer gas disk due to shear motions 
from the resettling gas. 

The dynamical models yield the 3D velocity vector of 
the galaxies, the peak ram pressures, and the times to 
peak ram pressure. In the case of a smooth, static, and 
spherical intracluster medium, peak ram pressure occurs 
during the galaxy's closest approach to the cluster center, 
i.e. when the galaxy's velocity vector is perpendicular to 
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its distance vector. Under these conditions the galaxy's 
present line-of-sight distance and its 3D position during 
peak ram pressure can be calculated. The linear orbital 
segments derived in this way together with the intracluster 
medium density distribution from Schindler et al. (1999) 
are consistent within a factor of 2 with the dynamical 
simulations for NGC 4501, NGC 4330, and NGC 4569. To 
resolve the discrepancy either a 2 times higher intraclus- 
ter medium density than derived from X-ray observations 
by Schindler et al. (1999) and/or a 2 times higher strip- 
ping efficiency than assumed by the dynamical models, is 
needed. The impact parameters of the galaxy orbits vary 
between 200 kpc and 600 kpc. Compared to NGC 4501, 
NGC 4330, and NGC 4569, NGC 4388 requires a still 2 
times higher local intracluster medium density or a slightly 
different direction of the galaxy's motion as derived from 
the dynamical model. 

It is suggested that NGC 4438 has moved within M 86's 
stripped hot X-ray gas during at least ~ 100 Myr, flew 
through the supposedly ram pressure stripped X-ray tail 
of M 86, recently left M 86's Mach cone, and is now en- 
countering the Virgo intracluster medium. In this scenario 
the trajectory of NGC 4438 is located ~ 200 kpc behind 
M 86 and the interaction with NGC 4435 is responsible 
for NGC 4438 's prominent tidal arms. 
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Appendix A: The region around M 86 

As shown in Sect. [H the proposed impact parameters 
of the two galaxies closest to M 86, NGC 4438 and 
NGC 4388, are not consistent with the peak ram pressure 
deduced from the dynamical models. These discrepancies 
can be due to the parameters deduced from the dynamical 
models. For NGC 4388 the local intracluster medium den- 
sity can well be enhanced by a factor of 2. I will explore 
the inconsistencies further in this Appendix. 
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A.l. The tidal interaction of NGC 4438 

The case of NGC 4438 is more complicated than that of 
NGC 4388. Dynamical models suggest that NGC 4438 had 
a gravitational interaction with NGC 4435 (Combes et 
al. 1988, VoUmer et al. 2005). Recent deep Ha observa- 
tions of the M 86 region (Kenney et al. 2008) revealed 
a highly complex and disturbed interstellar/intracluster 
medium. NGC 4438 is connected to M 86 by several faint 
Ha filaments, which might suggest a tidal interaction be- 
tween the two galaxies, as advocated by Kenney et al. 
(2008). The timescale and the relative galaxy velocities 
of the two scenarios (encounter with M 86 or NGC 4435) 
are about the same. The higher mass of M 86 compared 
to that of NGC 4435 allows a larger impact parameter 
for an NGC 4438 - M 86 encounter. Since M 86 has its 
own hot X-ray gas halo, the need for strong ongoing ram 
pressure that is presently affecting the ISM of the inner 
galactic disk, as shown in VoUmer et al. (2005, 2008), im- 
plies that ram pressure stripping due to M 86 's hot X-ray 
gas certainly did not remove, and might not even have 
significantly affected the ISM of NGC 4438. The recent 
discovery of a ~ 190 Myr-old starburst in NGC 4435 by 
Panuzzo et al. (2007) favors the interaction scenario be- 
tween NGC 4438 and NGC 4435. 

A. 2. The trajectory of NGC 4438 

In the following I assume that the faint Ha filaments 
connecting NGC 4438 to M 86 indicate NGC 4438 's 
direction of motion with respect to M 86. The 
model velocity vector of NGC 4438 (in kms^^) is 
= (-1440, -260, -1360il (Table [S]). According to 
Finoguenov et al. (2004) and RandaU et al. (2008), M 86 
has a significant velocity component towards the south. 

The 3D velocity vector of M 86 is vuse = 
{vx.Vy, -1500) kms~^ NGC 4438's velocity vector within 
the infalling subcluster of M 86 is thus At> — f N4438 — 
vmm = (-1440 -w^,0.2(-1440-Vx), 140) km^^, where 
the y-component is calculated to fit the direction of the 
faint Ha filaments connecting NGC 4438 to M 86. In 
this way I obtain Vy = -260 - 0.2(1440 -f v^) km~^ 
Following the Randall et al. model leads to Vx = Vy = 
-457 kms~^ M 86's total velocity with respect to M 87 
is thus 1630 kms~^. NGC 4438's relative velocity vector 
with respect to M 86 is Av = (-983,197,140) kms'^ 
Consequently, NGC 4438 has evolved until recently on a 
west-east trajectory within the hot X-ray gas of M 86 at 
a speed of ~ 1000 kms~^. 

Following the Mach cone geometry of Randall et al. 
(2008) and assuming that NGC 4438 is just leaving 
the Mach cone yields a line-of-sight distance between 
NGC 4438 and M 86 of - 200 kpc. This is smaller than 
the physical length of M 86's X-ray tail (380 kpc; Randall 
et al. 2008). 



^ negative numbers indicate directions to the left, bottom, 
and in front of M 87 



A. 3. Ram pressure stripping of NGC 4438 

To estimate the ram pressure exerted on NGC 4438 by 
M 86's hot X-ray gas, its gas density needs to be known. 
Randall et al. (2008) estimate the density in M86's bright 
X-ray plume to be n ~ 10"'^ cm~'^. The gas density at 
larger distances from M 86 might thus be a few 10~^ cm~^. 
On the other hand, the model deduced ram pressure for 
NGC 4438 (VoUmer et al. 2005) yields an intracluster 
medium density of ~ 10~^ cm~^ if a galaxy velocity of 
2000 kms~^ is assumed. 

In the proposed scenario NGC 4438 has been mov- 
ing during the last 100 Myr through the hot X-ray gas of 
M 86 (n ~ 10^^ cm^'^), has just recently (a few 10 Myr 
ago) left M86's Mach cone (see Fig. 13 of Randall et al. 
2008 for an illustration), and is presently encountering 
the Virgo cluster intracluster medium {n ^ 10^"^ cm"'^). 
Consequently, NGC 4438 is located at the same line-of- 
sight distance as M 86, i.e. 1.0 ± 0.8 Mpc behind M 87 
(Mei et al. 2007). At the position of NGC 4438 Bohringer 
et al. (1994) have observed a filament of intracluster gas. 
Finoguenov et al. (2004) estimate the gas density in the 
filament to be ~ 10"'^ cm"-^. 

I conclude that NGC 4438 is not stripped by the spher- 
ical component of the Virgo intracluster medium. The 
stripping rather occurs at a physical distance from M87 
between 500 kpc, assuming a line-of-sight distance be- 
tween M 86 and M 87 of 200 kpc (RandaU et al. 2008), and 
1 Mpc, adopting the M 86-M 87 distance from Mei et al. 
(2007). Together with NGC 4522, NGC 4438 is the second 
galaxy which is undergoing strong ram pressure stripping 
at a location in the Virgo cluster where this would not 
be expected from a classical model assuming a spherical, 
smooth, and static intracluster medium distribution. 

A.4. The ram pressure stripped gas tail of NGC 4388 

Given that the 100 kpc long Hi taU of NGC 4388 is ex- 
ceptional^ for a galaxy at a projected distance from M 87 
less than 0.5 Mpc, I also consider a possible ram pressure 
stripping event with M 86's intracluster medium within 
the Mach cone ^ 200 Myr ago as proposed by Oostcrloo & 
van Gorkom (2005). This would explain why such long Hi 
tails are not observed in the other galaxies of the present 
sample. An argument against this hypothesis is the ex- 
treme difference between the radial velocities of NGC 4388 
and M 86 (2700 kms~^) compared to the relatively low 
tangential velocity (80 kpc/150 Myr=500 kms~^) which 
implies a ratio between the line-of-sight and the projected 
tail length of ^ 5 and thus a physical length of ^ 400 kpc. 
Another possibility might be that the gravitational po- 
tential of M 86 enhances locally the intracluster medium 
density near M 87 and the embedded gas of NGC 4388's 
stripped tail, leading to an enhanced recombination of the 
ionized gas in the tail which makes it observable in Hi. 

There are two other prominent gas tails detected in 
NGC 4532/DDO 137 (Koopmann et al. 2008) and NGC 4254 
(Haynes et al. 2007). Both are most probably of tidal nature. 



